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The  turbulence  structure  of  convective  heat  transfer  was  studied  experimentally  in  complex  three-  ; 
dimensional  and  separated  turbulent  boundary  layers.  Three  test  cases  whose  fluid  dynamics'  have  been 

well  documented  were  examined.  In  case  1,  time-  resolved  surface  heat  transfer  was  measured  in  the  nose 
region  of  a  wing-body  junction  formed  by  a  wing  and  a  flat  plate.  Mean,  statistical,  and  spectral 
characteristics  of  the  surface  heat  transfer  are  reported.  The  effects  of  wing  shape  were  investigated  by 
measuring  the  surface  heat  transfer  in  the  nose  region  of  a  modified  NACA  0020,  a  streamlined  cylinder 
shape  and  NACA  0015.  The  effectiveness  of  a  flow  control  device  to  reduce  surface  heat  transfer  is 
reported.  In  case  2,  simultaneous  surface  heat  flux  and  temperature  protiles  were  measured  at  8  locations 
in  the  spatially-developing  pressure-driven  three-dimensional  turbulent  boundary  layer  upstream  of  the 
wing-body  junction.  In  case  a,  simultaneous  surface  heat  flux  and  temperature  profiles  were  measured  at 
18  stream-wise  locations  in  a  mean  2-dimensional  adverse-pressure-gradient  separating  turbulent  boundary 
layer.  Mean,  statistical  and  spectral  heat  flux  and  temperature  data  are  reported.  Mean  ejection  frequencies, 
turbulence  length  scales,  inclination  angles  of  the  turbulence  structure,  and  coherency  between  the  inner 
and  outer  regions  of  the  flow  were  examined  from  these  results.  Several  useful  correlations  between 
surface  heat  transfer  and  velocity  are  reported. 
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1.  Objectives  of  the  Research  Effort 

It  is  well  known  that  in  many  turbulent  separated  flow  cases  the  convective  heat 
transfer  rate  is  much  greater  than  for  an  attached  flow.  While  there  are  global  heat 
transfer  measurements  for  such  cases,  no  detailed  study  of  the  relationships  between 
the  fluid  turbulence  and  the  heat  transfer  have  been  made. 

Here  simultaneous  time-  and  spatially-resolved  heat  transfer  and  temperature 
profile  measurements  were  made  in  a  program  that  began  6/16/91  for  three  low  speed 
separated  flow  test  cases  for  which  detailed  laser  velocimeter  flowfield  measurements 
were  already  available: 

(1)  The  3-D  steady  free-stream  velocity  wing/body  junction  flow  of 
Devenport  and  Simpson  (1987,  1990),  which  demonstrates  a  self-induced 
bimodal  unsteady  separated  flow  at  the  nose  that  undoubtedly  is 
responsible  for  the  high  heat  transfer  rates  observed  there. 

(2)  The  pressure-driven  3-D  turbulent  boundary  layer  flow  of  Okjmen  and 
Simpson  (1994)  which  exhibits  an  anisotropic  eddy  viscosity  and 
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reduction  in  coherence  between  the  inner  and  outer  regions  due  to  three- 
dimensionality. 

(3)  The  mean  2-D  steady  free-stream  velocity  adverse-pressure-gradient- 
induced  separated  flow  of  Simpson,  Chew,  and  Shivaprasad  (1981  a,  b), 
which  was  a  computational  test  case  for  the  1980-81  APOSR-Stanford 
Conferences  on  Complex  Turbulent  Flows. 

For  the  first  case,  several  passive  controls  of  the  bimodal  unsteady  flow  were  also  examined  in 
regard  to  their  effectiveness  in  inhibiting  heat  transfer. 

The  results  of  the  research  on  all  three  cases  is  of  value  to  both  basic  and  applied  issues. 
The  basic  issues  of  turbulent  transport  of  a  passive  scaler  (temperature)  in  separated  flow  is  of 
primary  scientific  merit.  The  first  case  is  of  interest  because  of  high  heat  transfer  at  the  blade/hub 
junction  of  jet  engine  turbines.  The  last  two  cases  have  airfoil  and  axial  compressor  types  of 
pressure  distributions. 


2.  Description  of  Test  Cases 

2.1  Case  1  -  Wing/Body  Junction  Flow  Case 

Figure  1  shows  a  schematic  of  the  3:2  elliptic  nose,  NACA  0020  airfoil  tail  wing  has  been 
used  in  experimental  studies  of  the  flowfield  and  turbulence  for  this  type  of  pressure-gradient- 
driven  3-D  flow  (Devenport  and  Simpson,  1987,  1988,  1990).  (This  flow  is  being  used  as  a  test 
case  for  the  1990-91  Collaborative  Testing  of  Turbulence  Models.)  When  the  approach  boundary 
layer  encounters  the  wing,  it  separates  in  front  of  the  nose.  Closer  to  the  nose  a  "horseshoe" 
vortex  forms  which  is  stretched  around  the  body.  The  separation  line  (converging  surface  skin 
friction  lines)  stretches  around  the  sides  of  the  wing. 


Devenport  and  Simpson  observed  bimodal  velocity  probability  histograms  in  the  zone 
outlined  by  the  solid  line.  Further  research  showed  that  the  flow  in  this  zone  aperiodically 
switched  from  one  mode  to  another.  Turbulence  energy  production  and  turbulent  stresses  an 
order  of  magnitude  higher  than  upstream  were  found  in  this  zone  because  of  this  low  frequency 
self-induced  switching. 


Hydrogen-bubble  flow  visualization  experiments  have  been  done  on  a  2  foot  chord  clear 
plexiglas  model  in  a  water  tunnel.  Figure  2  shows  a  schematic  of  the  self-induced  unsteady 
phenomena  revealed  from  high  speed  videos. 

At  time  t,  a  large  horseshoe  vortex  exists  in  front  of  the  nose,  produced  by  high  velocity 
free-stream  fluid  impinging  on  and  moving  down  the  nose.  Because  the  vortex  lines  of  this  flow 
are  stretched  around  the  wing,  the  cross-sectional  area  of  the  vortex  decreases  at  increasing  times 
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t,,  t3,  and  t4.  Meanwhile,  a  separation  vortex  forms  downstream  of  separation,  increasing  with 
circulation  strength  at  increasing  times  t3  and  t4.  At  some  time,  the  separation  vortex  either 
merges  with  the  front  of  the  horseshoe  vortex  or  moves  up  over  the  horseshoe  vortex  in  "leap 
frog"  fashion  before  merging  with  the  horseshoe  vortex.  The  resulting  merger  creates  a  new 
large  horseshoe  vortex  and  the  aperiodic  process  begins  anew. 

Each  repetition  of  this  sequence  is  not  exactly  repeatable  either  in  spatial  location  or 
period  of  events.  However,  we  believe  that  this  is  a  vorticity-dominated  flowfield  that  is 
controlled  by  the  stretching  of  the  horseshoe  vortex  and  the  separation  vortex.  Thomas  (1987) 
and  Cantwell  (1990)  observed  self-induced  unsteadiness  and  merging  of  horseshoe  and  separation 
vortices  in  two  independent  experiments  of  a  laminar  boundary  layer  approaching  a  circular 
cylinder. 

Furthermore,  we  believe  that  this  bimodal  aperiodic  phenomenon  accompanies  all 
wing/body  junction  flows  and  is  responsible  for  observed  high  turbulence  intensities,  high 
pressure  fluctuations,  and  high  heat  transfer  rates.  This  is  of  great  importance  to  gas  turbine  heat 
transfer,  where  hub  and  casing  boundary  layers  produce  large  heat  transfer  rates  at  the  junctions 
with  blades. 

2.2.  Case  2  -  Pressure-Driven  3-D  Turbulent  Boundary  Laver  Test  Case 

A  three-dimensional  boundary  layer  is  created  when  a  spanwise  pressure  gradient  or  wall 
shear  stress  cause  the  flow  direction  to  vary  across  the  boundary  layer  as  shown  in  Fig.  3.  Three- 
dimensional  turbulent  boundary  layers  occur  in  many  practical  applications  where  convective  heat 
transfer  is  significant.  Flow  in  gas  turbines  is  one  well-known  example.  Previous  research  has 
shown  that  three-dimensionality  reduces  the  turbulent  transport  of  momentum  and  heat  and  that 
the  eddy  viscosity  is  not  isotropic. 

The  test  case  used  here,  shown  in  Fig.  4,  differs  from  many  other  experimentally 
examined  three-dimensional  flows  in  that  the  mean  flow  variables  depend  on  three  spatial  axes 
rather  than  two  axes,  such  as  flows  in  which  the  three-dimensionality  of  the  flow  has  been 
generated  either  by  a  rotating  cylinder  or  by  a  pressure  gradient  in  one  direction  only  throughout 
the  flow.  Initially,  the  spanwise  pressure  gradient  created  by  the  presence  of  the  wing  turns  the 
flow  away  from  the  wing.  Further  downstream,  the  spanwise  pressure  gradient  changes  sign  and 
the  flow  is  turned  back  toward  the  wing.  Both  adverse  and  favorable  stream-wise  pressure 
gradients  also  occur.  The  turbulence  structure  is  strongly  three-dimensional.  This  well- 
documented  flowfield  can  be  used  as  a  computational  test  case  of  turbulence  models  since  the 
initially  mean  2-D  flow  upstream  develops  into  a  non-equilibrium  3-D  flow  and  then  relaxes  back 
toward  a  mean  2-D  flow  far  downstream. 

2.3.  Case  3  -  Mean  2-D  Steady  Free-Stream  Velocity  Adverse-Pressure-Gradient-Induced 
Separated  Flow 

The  airfoil-type  stream-wise  free-stream  velocity  flow  of  Simpson,  Chew,  and  Shivaprasad 
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(1981  a,  b)  on  a  flat  wall  (Figure  5)  was  a  computational  test  case  for  the  1980-81  AFOSR- 
Stanford  Conferences  on  Complex  Turbulent  Flows  and  the  1990-91  Collaborative  Testing  of 
Turbulence  Models.Numerous  Codes  and  turbulence  models  have  attempted  to  calculate  this  flow, 
resulting  in  some  improved  modelling  (simpson,  1985,  1986).  The  emphasis  in  the  author’s  group 
has  been  in  determining  the  flow  physics  from  the  experimental  data. 

Figure  6  briefly  summarizes  just  a  few  of  the  distinct  features  of  the  type  of  separated 
flow.  Turbulent  diffusion  of  turbulence  kinetic  energy  controls  the  backflow  rate  (Agarwal  and 
Simpson,  1989).  Diffusion  and  dissipation  are  dominant  processes  in  the  backflow  while  the 
Reynolds  shear  stress  -puv  and  turbulence  energy  production  are  small.  The  mean  backflow 
velocity  profile  of  Simpson  (1983)  holds. 

3.  Experimental  Facilities  and  Instrumentation 

3.1.  Boundary  Laver  Tunnel 

The  wind  tunnel  is  an  open-circuit  type  and  is  powered  by  a  centrifugal  blower.  Air  from 
the  blower  is  supplied  to  the  test  section  after  passing  through  a  fixed-setting  damper,  a  plenum, 
a  section  of  honeycomb,  7  screens  which  are  used  to  remove  much  of  the  turbulence  intensity, 
and  a  4: 1  contraction  nozzle  to  further  reduce  the  turbulence  and  to  accelerate  the  flow  to  test 
speed.  The  potential  core  of  the  flow  entering  the  test  section  is  uniform  to  within  0.5%  in  the 
span  wise  direction  and  1%  in  the  vertical  direction  with  a  turbulence  intensity  of  0.1%  at  27  m/s.1 

3.1.2.  Wing/Body  junction  Flow 

Figure  7  is  a  side  view  of  the  6-m-long  and  0.91-m-wide  test  section.  The  upper  wall  is 
made  from  Plexiglas®  reinforced  with  aluminum  channel.  The  glass  side  walls  are  lined  internally 
with  removable  6.4  mm  thick  Plexiglas  sheets.  The  lower  wall  is  made  from  19  mm  thick  fin- 
form  plywood.  Flow  entering  the  test  section  is  subjected  to  a  further  1.5:1  contraction  produced 
by  the  shape  of  the  upper  wall.  A  throat  is  reached  1.63  m  downstream  of  the  entrance  where 
the  section  is  254  mm  in  height.  Downstream  of  the  throat  the  upper  wall  is  almost  parallel  to 
the  flat  lower  wall,  diverging  gradually  from  it  with  distance  downstream  to  account  for  boundary 
layer  growth. 

In  order  to  create  a  new  fluid  boundary  layer  at  the  beginning  of  the  heat  flux  surface, 
a  heated  false  floor  was  placed  in  the  tunnel  and  a  suction  slot  was  opened  in  the  tunnel  floor 
as  shown  in  Fig.  7.  The  upstream  boundary  layer  formed  on  the  tunnel  floor  is  sucked  out 
through  the  slot  and  a  new  boundary  layer  develops  from  the  rounded  leading  edge  of  the  false 
floor.  The  thin  boundary  layer  is  tripped  by  a  12.7  mm  wide  strip  of  120-grade  sandpaper, 
located  with  its  leading  edge  70  mm  downstream  of  the  false  floor  leading  edge.  The  forward 
edge  of  the  suction  slot  is  located  145  mm  downstream  of  the  throat  of  the  test  section.  The  false 
floor  rests  with  its  upper  surface  35  mm  above  the  tunnel  floor  and  its  leading  edge  116  mm 
downstream  of  the  throat.  The  false  floor  was  made  of  3  sections  of  1.59  cm  thick  aluminum 
plate.  A  trailing  edge  flap  set  at  an  angle  of  7°  was  used  to  produce  a  constant  static  pressure  on 
the  tunnel  floor  throughout  the  heated  test  section. 
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The  wing  was  mounted  in  the  test  section  at  zero  incidence  and  sweep  with  its  leading 
edge  1.17m  downstream  of  the  leading  edge  of  the  false  floor.  Since  the  measurements  were 
made  only  with  the  wing  at  a  zero  angle  of  attack,  it  was  possible  to  reduce  blockage-induced 
pressure  gradients  by  contouring  the  side  walls  of  the  wind  tunnel  to  approximately  follow  the 
streamlines  produced  by  the  3:2  elliptic  nose/NACA  0020  tailed  wing  in  unbounded  potential 
flow.  This  was  accomplished  by  removing  the  6.4  mm  thick  Plexiglas  sheets  which  line  the 
tunnel  sidewalls  from  the  region  surrounding  the  wing.  This  effectively  increased  the  width  of 
the  test  section  by  12.7  mm  from  a  location  330  mm  upstream  of  the  wing  leading  edge  to 
another  203  mm  downstream  of  its  trailing  edge.  The  abrupt  comers  at  the  edges  of  the  liner 
were  faired  over  using  adhesive  tape.  This  same  arrangement  was  used  for  all  wing  shapes. 

The  aluminum  false  floor  was  heated  and  held  at  a  constant  and  uniform  temperature. 
Silicone-rubber-insulated  electric  resistance  heaters  were  held  against  the  bottom  of  the  false  floor 
by  1.9  cm  thick  styrofoam  sheet  which  was  secured  by  bolts.  T-type  thermocouples  were  soldered 
into  3.18  mm  diameter  brass  tube  and  press  fit  into  the  aluminum  floor  for  feedback  control. 
Power  to  the  heaters  was  controlled  by  5,  Eurotherm  810,  3-mode  process  controllers,  each 
connected  to  a  Eurotherm  831  SCR  power  supply.  Because  of  the  high  thermal  conductivity  of 
aluminum  and  the  thickness  of  the  false  floor,  the  surface  was  nearly  isothermal.  The  automatic 
controllers  held  the  surface  temperature  constant  and  uniform  to  within  ±0.5°C.  The  virtual  origin 
of  the  turbulent  boundary  layer  was  located  1  cm  downstream  of  the  leading  edge  of  the  false 
floor.  The  unheated  starting  length  was  approximately  5  cm  and  was  therefore  negligible. 

The  heat  flux  probe  was  mounted  in  a  cam  system  installed  in  the  false  floor  in  front  of 
the  wing.  The  cams  were  made  of  1.59  cm  thick  aluminum  plate  and  were  heated  from  beneath 
by  electric  resistance  heaters.  Rotating  one  or  both  cams  allowed  the  probe  to  be  positioned 
throughout  a  12.7  cm  radius  region  in  front  of  the  wing.  The  cam  system  is  illustrated  in  Figure 
7  of  Lewis  et  al.  (1993). 

3.1.2.  Pressure-Driven  3-D  Turbulent  Boundary  Laver  Boundary  Laver 

To  increase  the  momentum  Reynolds  number  to  match  the  flow  conditions  of  Oilmen  and 
Simpson  (1995),  the  heated  false  floor  used  in  the  wing/body  junction  case  was  replaced  with  the 
permanent  heatable  floor  shown  in  Figure  8.  The  floor  of  the  tunnel  was  heated  from  a  distance 
121.9  cm  downstream  of  the  test  section  entrance  (41.07  cm  upstream  of  the  throat)  for  a 
distance  of  335.3  cm.  The  heated  floor  is  made  from  six  sections.  The  middle  three  sections 
consist  of  the  floor  used  in  the  wing/body  junction  case.  The  three  new  sections  are  made  from 
0.375  thick  aluminum  plate.  As  before,  the  floor  is  heated  from  beneath  with  silicone  rubber 
heaters.  To  ensure  better  thermal  contact  between  the  heaters  and  the  aluminum  plates,  the 
heaters  were  epoxied  to  the  bottom  of  the  floor  using  Silicone  RTV  epoxy.  The  cam  system  was 
replaced  with  a  heated  plate  with  holes  milled  at  the  eight  locations  where  Oilmen  and  Simpson 
(1995)  made  LDV  measurements.  The  heat  flux  gage  is  press-fit  into  these  holes.  A  ramp  was 
installed  upstream  of  the  aluminum  floor  to  gradually  adjust  from  the  height  of  the  leading  edge 
of  the  test  section  entrance  to  the  height  of  the  raised  floor.  The  ramp  is  121.9  cm  long  and  3.493 
cm  high.  The  ramp  consists  of  a  plexiglass  surface  supported  by  an  aluminum  frame. 

The  upper  wall  of  the  test  section  was  adjusted  to  give  the  same  stream-wise  velocity 
variation  as  Oilmen  and  Simpson  (1995).  The  corner  gaps  between  the  upper  wall  and  glass  side 
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walls  were  covered  with  a  flexible  polyurethane  plastic  sheet  to  prevent  leakage  at  these  comers. 

3.1,3.  Mean  2-D  Steady  Free-Stream  Velocity  Adverse-Pressure-Gradient-Induced  Separated  Flow 

For  this  test  case,  the  inner  liners  of  the  test  section  was  removed  and  the  upper  wall  was 
adjusted  as  shown  in  Fig.  4.  An  active  boundary  layer  control  system,  described  by  Simpson, 
Chew,  and  Shivaprasad  (1980),  is  used  to  eliminate  preferential  separation  of  the  curved-top-wall 
boundary  layer.  Highly  two-dimensional  wall  jets  of  high-velocity  air  are  introduced  at  the 
beginning  of  each  of  the  244-cm-long  sections.  At  the  latter  two  stream-wise  locations  the 
oncoming  boundary  layer  is  partially  removed  by  highly  two-dimensional  suction  systems. 

3.2  Heat  Flux  Microsensor. 


Heat  flux  measurements  were  made  using  both  a  high-  frequency-response,  low  output  gage  and 
a  high  sensitivity  gage.  The  high-frequency-response  gage  was  a  Heat  Flux  Microsensor 
manufactured  by  Vatell  Corp.  Similar  gages  have  been  used  in  several  other  studies. of  turbulent 
heat  flux.  The  Heat  Flux  Microsensor  is  a  layered  type  gage  which  operates  by  relating  the  heat 
flux  to  the  measured  temperature  drop  across  a  thin  thermal  resistance  layer.  The  Microsensor 
has  a  flat  frequency  response  up  to  50  kHz  and  outputs  a  voltage  directly  proportional  to  the  heat 
flux.  The  thin  construction  of  the  gage  (<2pm)  is  physically  unobtrusive  in  the  flow  and  causes 
minimal  disruption  of  the  thermal  boundary  layer  due  to  a  step  change  in  temperature.  The 
sensing  area  of  the  gage  used  in  the  wing/body  junction  case  was  3  mm  by  3  mm.  The  sensing 
area  of  the  gage  is  used  in  case  2  and  case  3  was  4  mm  by  6.5  mm.  The  gage  was  always 
aligned  so  that  the  4-mm  length  was  in  the  stream-wise  direction.  The  output  voltage  of  the 
Microsensor  was  amplified  1000  times  by  a  differential  dc-powered  amplifier  designed  and  built 
by  Vatell  Corp. 

3.3.  Schmidt-Boelter  gage. 

Due  to  the  low  output  voltage  of  the  Microsensor  used  in  case  1,  a  second,  high 
sensitivity  Schmidt-Boelter  gage  manufactured  by  Medtherm  Corp.  (Model  #  8-2-. 625-36- 
20893T)  was  used  to  measure  the  mean  time-averaged  heat  flux.  The  Schmidt-Boelter  gage  has 
an  output  voltage  directly  proportional  to  the  heat  flux  with  a  sensitivity  of  3.9  ±  0.2 
mV/(W/cm2)  and  a  response  time  of  approximately  250  msec.  The  diameter  of  the  sensing  area 
of  the  gage  is  3.18  mm.  Temperature  of  the  gage  was  measured  by  a  T-type  thermocouple 
mounted  inside  the  gage  by  the  manufacturer.  The  output  voltage  from  the  Schmidt-Boelter  gage 
was  amplified  1000  times  filtered  by  a  300  Hz  low-pass  filter. 

3.4.  Cold  Wire 


Temperature  was  measured  using  a  hot-wire  anemometer  operating  in  the  constant-current  mode. 
A  Dantec  56C01  anemometer  with  a  Dantec  56C20  constant-temperature  bridge  was  used.  The 
probe  used  is  a  TSI  1261 A  -  P.5  miniature  boundary  layer  probe  which  has  been  modified  by 
TSI.  The  original  TSI  1261 A  -  T1.5  probe  was  modified  by  replacing  the  5  pm  diameter  x  1.5 
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mm  length  tungsten  wire  with  a  1 .25  pm  dia.  x  0.5  mm  length  platinum  wire.  The  small  diameter 
wire  provides  greater  frequency  response  due  to  its  small  thermal  inertia.  The  smaller  length  of 
the  wire  reduces  spatial  averaging  of  small  turbulent  structures.  The  high  conductivity  of  platinum 
increases  the  frequency  response  of  the  wire  as  well  as  enabling  the  length  of  the  wire  to  be 
reduced  without  increasing  heat  loss  to  the  support  needles.  The  frequency  response  of  the  wire 
was  measured  by  both  an  external  heating  method  and  a  current  injection  method  (Lewis,  1995). 
The  wire  frequency  recency  was  flat  up  to  approximately  2  kHz  at  all  velocities  and  turbulence 
intensities  of  interest. 

3.5.  Data  Acquisition 

For  case  1,  The  heat  flux  signal  was  sampled  by  an  HP  3562A  Dynamic  Signal  Analyzer 
at  a  frequency  of  10240  hz.  One  hundred  blocks  of  2048  samples  were  taken  over  a  period  of 
approximately  1  min  and  recorded  on  floppy  disk.  Digital  data  reduction  was  performed  later. 
Time-  and  spatially-resolved  surface  heat  flux  measurements  were  made  in  the  nose  region  of  a 
wing/body  junction.  Both  mean  and  fluctuating  heat  fluxes  were  measured  at  60  locations  in  the 
nose  region  of  the  3:2  elliptic  nose/NACA  0020  tail  airfoil.  Measurements  were  also  made  along 
the  line  of  symmetry  in  front  of  two  other  wing  shapes  to  examine  the  effects  of  geometry. 

In  addition,  since  we  have  examined  several  passive  control  devices  to  eliminate  or  reduce 
the  bimodal  flow  switching,  we  examined  the  benefit  of  one  of  these  passive  control  devices  in 
reducing  heat  transfer  around  the  junction  for  three  different  wing  shapes. 

For  case  2,  simultaneous  temperature  and  surface  heat  flux  signals  were  taken  by  the  HP 

3564  A.  60  loads  of  2048  samples  of  each  signal  were  taken  at  a  frequency  of  10240  hz  and 

recorded  on  floppy  disk.  Spectral  analysis  was  performed  in  real  time  by  the  HP  3562A.  200 
ensemble  averages  of  the  heat  flux  and  temperature  power  spectra,  co-spectra,  and  spectral 
coherence  were  calculated  at  a  sampling  frequency  of  10240  hz  and  recorded.  Simultaneous 
measurement  of  temperature  profiles  and  surface  heat  flux  were  made  at  8  locations  in  the  three- 
dimensional  boundary  layer  upstream  of  the  wing/body  junction  outside  of  the  region  near  the 
wing  where  large  mean  flow  stream-wise  vortices  are  dominant.  At  each  profile  location,  surface 
heat  transfer  was  recorded  simultaneously  with  temperature  at  30  locations  spaced  logarithmically 
across  the  boundary  layer. 

For  case  3,  simultaneous  temperature  and  surface  heat  flux  signals  were  taken  by  the  HP 

3564  A.  60  loads  of  2048  samples  of  each  signal  were  taken  at  a  frequency  of  409.6  hz  and 

recorded  in  real  time  on  floppy  disk.  Spectral  analysis  was  performed  in  real  time  by  the  HP 
3562A.  1500  ensemble  averages  of  the  heat  flux  and  temperature  power  spectra,  co-spectra,  and 
spectral  coherence  calculated  at  a  sampling  frequency  of  10240  hz  were  recorded.  Simultaneous 
measurement  of  time-resolved  surface  heat  transfer  and  temperature  were  measured  at  30 
logarithmically  spaced  y-locations  at  each  of  14  stream-wise  locations  from  the  beginning  of  the 
adverse  pressure  gradient  through  separation. 
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4.  Catalog  of  Experimental  Results 


4.1.  Wing/Bodv  Junction 

1)  Mean  surface  heat  transfer 

2)  Variance, skewness  and  flatness  of  surface  heat  transfer  fluctuations 

3)  probability  density  functions  of  surface  heat  transfer  fluctuations 

4)  spectra  of  surface  heat  transfer  fluctuations 

5)  St'  correlation  relating  mean  surface  heat  transfer  to  turbulence  intensity 

6)  Effects  of  nose  shape  on  mean  surface  heat  transfer 

7)  Effects  of  flow  control  fence  on  mean  and  variance  of  surface  heat  transfer 

4.2.  Pressure-Driven  3-Dimensional  Turbulent  Boundary  Layer 

1)  Profiles  of  mean  temperature 

2)  mean  surface  heat  transfer 

3)  profiles  of  variance,  skewness,  and  flatness  of  temperature  variance 

4)  spectra,  cospectra  and  coherence  of  temperature  and  surface  heat  transfer  fluctuations. 

5)  Coherence  length  scales 

6)  Estimates  of  mean  angles  of  large-scale  turbulent  structures 

7)  Correlation  of  velocity  and  temperature  spectra. 

8)  Mean  ejection  frequencies. 

9)  Beta  correlation  relating  temperature  fluctuations  to  velocity  fluctuations. 

10)  Abrahamson  and  Eaton  Correlation  relating  mean  surface  heat  transfer  to  magnitude 
of  enthalpy  thickness  vector. 

4.3.  Mean  2-D  Steady  Free-Stream  Adverse-Pressure-Gradient-Induced  Separated  Flow 

1)  Profiles  of  mean  temperature 

2)  mean  surface  heat  transfer 

3)  profiles  of  variance,  skewness,  and  flatness  of  temperature  variance 

4)  spectra,  cospectra  and  coherence  of  temperature  and  surface  heat  transfer  fluctuations. 

5)  Coherence  length  scales 

6)  Estimates  of  mean  angles  of  large-scale  turbulent  structures 

7)  Correlation  of  velocity  and  temperature  spectra. 

8)  Mean  ejection  frequencies. 

9)  Beta  correlation  relating  temperature  fluctuations  to  velocity  fluctuations. 

10)  Abrahamson  and  Eaton  Correlation  relating  mean  surface  heat  transfer  to  magnitude 
of  enthalpy  thickness  vector. 

11)  St' correlation 
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5.  Summary  of  Unique  Experimental  Results 


This  study  is  unique  among  heat  transfer  studies.  To  the  authors’  knowledge,  this  is  the 
only  study  of  statistical  and  spectral  properties  of  temperature  and  surface  heat  transfer  in  three- 
dimensional  and  separated  turbulent  flows.  A  unique  feature  of  the  wing/body  junction  flow  case 
is  the  use  of  several  different  wing  shapes  to  evaluate  the  effects  of  wing  shape  on  surface  heat 
transfer.  The  second  case  is  the  only  known  investigation  of  the  turbulence  structure  of 
convection  in  a  three-dimensional  turbulent  boundary  layer.  The  test  cases  studied  here  have  well- 
documented  fluid  dynamics  and  are  suitable  test  cases  for  validation  of  CFD  codes  and 
turbulence  models. 

6.  Future  Work 


The  dissertation  will  be  completed  in  late  summer.  Two  journal  articles  are  being  prepared 
from  the  results  from  the  3-D  case  and  the  2-D  separating  case. 
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Figure  1 


Sketch  of  the  wing-body  junction  showing  the  horseshoe  vortex. 


1.  HORSESHOE  VORTEX  EOFWS  AT  THE  NOSE  OF  THE  WING  AS  SHOWN  AT  TIME  T^. 

2.  AREA  OF  HORSESHOE  VORTEX  DECREASES  AS  VORTEX  IS  STRETCHED  AROUND  THE 
NOSE  OF  THE  OF  THE  WING  AT  TIMES  T2  AND  LATER. 

3.  SEPARATION  VORTEX  FOIWS  AND  GROWS  AT  TIMES  T_  AND  T„. 

3  4 

4.  AT  SOME  TIME,  THE  SEPARATION  VORTEX  EITHER  MERGES  WITH  THE 
HORSESHOE  VORTEX  FROM  THE  FRONT  OR  MOVES  PARTLY  UP  AND  OVER  THE 
HORSESHOE  VORTEX  IN  A  LEAPFROG  FASHION  BEFORE  MERGING. 

5.  SE3JENCE  OF  EVENTS  REPEATS  IN  AN  APERIODIC  MANNER. 


Figure  2  Schematic  of  the  sequence  of  pairing  of  horseshoe  and  separation  vortices. 
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Figure  3  Sketch  of  typical  three-dimensional  boundary  layer  profile. 
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Figure  5  Sideview  schematic  of  the  test  section  with  the  steady  free-stream  separating 
turbulent  boundary  layer  on  the  bottom  wall. 


15 


A  FLOW  MODEL  WITH  THE  COHERENT  STRUCTURES  SUPPLYING 
THE  SMALL  MEAN  BACKFLOW.  ID  DENOTES  INCIPIENT  DETACHMENT; 
ITD  DENOTES  INTERMITTENT  TRANSITORY  DETACHMENT;  D  DENOTES 
DETACHMENT.  THE  DASHED  LINE  DENOTES  U  =  0  LOCATIONS. 

FROM  SIMPSON  FT  At  .  (1981b). 


1.  LARGE  EDDIES  GROW  DURING  DETACHMENT 

2.  URGE  EDDIES  SUPPLY  TURBULENCE  ENERGY  TO  BACKFLOW  AND  CONTROL 
OUTER  REGION  ENTRAINMENT  RATE. 

3.  URGE  EDDY  BEHAVIOR  SCALES  ON  MAXIMUM  SHEAR  STRESS. 

N.  DIFFUSION  AND  DISSIPATION  OF  TURBULENCE  ENERGY  IN  BACKFLOW. 

5.  SMALL  -uv  IN  BACKFLOW;  COLES  "UW-OF-THE-WALL"  DOES  NOT  APPLY. 


Figure  6  Instantaneous  backflow  behavior  structure  and  nature  of  backflow  important  for 
proper  modeling. 
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Upper  wall 


Figure  8  Side  view  of  the  boundary  layer  test  section  configuration  used  tor  the  pressure- 
driven  3-dimensional  turbulent  boundary  layer  test  case. 
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Time-  and  spatially-resolved  heat  tlux  measurements  are  reported  for  the  endwall  surface  in  the  turbulent. 

incompressible  flow  in  the  nose  region  of  a  wing-body  junction  formed  by  a  wing  ^cylindrical 

wing  and  the  flat  plate  were  heated  and  held  at  a  constant  and  uniform  temperature.  The  effects  cy  md 
IZ  tX on  heat  dux  were  investigated  by  taking  heat  dux  ™-uremen«  m .the 
eiliotic  nos«/NAC  A  0020  tail  shape,  a  circular  cylinder  with  a  wedge  tail,  and  an  t  AC.  • 
were  in^r^  upto  a  factor  of  3  over  ,h.  heat  (lux  rates  in  the  approach  boundary  layer.  The  rms  of  the  heat 
..  .  Tea  #  the  [nean  heat  flux  in  the  vortex-dominated  nose  region.  Away  from 

SxTng  i  hy  large-scale  unsteadiness  of  the  vortex.  A  new  three-dimensional  extensmn  *»«■»< 
correlation  is  proposed  to  account  for  the  effect,  of  the  horseshoe  vortex  on  heat  transfer  tn  the  regton.  Adjace 
to  the  wing  the  augmentation  in  heat  dux  is  due  to  a  change  in  the  mean  velocity  tie  d. 


Nomenclature 

C„  =  specific  heat 

£  =  output  voltage  at  gauge 

f  9  frequency,  Hz 

G(t(f)  -  auto  spectral  density  function  of  q 

h  =  convection  film  coefficient,  qi(T<  —  7'») 

lc  =  turbulent  kinetic  energy 

Pr  =  Prandtl  number 

q  -  time-mean  heac  flux 

q'  =  fiuctuating  component  of  heat  flux 

pe  =  Reynolds  number  based  on  (  ) 

5  '  =  heat  flux  gauge  sensitivity 

St  -  Stanton  number,  kt(pCpUrci) 

$('  =  Stanton  prime  correlation,  hiipC^u  rma.mux) 

T  -  temperature 

T,  =  film  temperature,  (T,  +  Tji2 

T.  =  gauge  surface  temperature 

7^  =  temperature  of  treestream  air 

r  *  =  maximum  thickness  of  a  given  wing 

U%  v,  W  -  time-averaged  velocity  components  in  x%  y,  r 
directions,  respectively 
(j  =  freestream  velocity 

Llrvv  w  =  fluctuating  velocity  components  in  .c,  y,  r 
directions,  respectively 

x,  v,  z  «  Cartesian  coordinate  system  axis 
xu'J  ~  =*  distance  measured  along  line  of  symmetry  of 

wing  from  wing  leading  edge  to  line  of  low 
shear 

x  =  distance  measured  along  line  of  symmetry  or 

wing  from  wing  leading  edge  to  line  of 
separation 

\T  =  temperature  difference,  7^  - 

=  gauge  emissivuy 
q  =  momentum  thickness 

p  =  density 
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=  Stefan  Boltzmann  constant. 
5.67  x  10*  W/(m-  FC) 


Subscripts 

rms 

rms, max 


root  mean  square  of  fluctuating  component 
maximum  rms  value  of  fluctuating  component 
present  in  the  boundary  layer  at  a  given  x.  : 
location 


Introduction 

4  WING/BODY  junction  flow  is  formed  when  the  bound- 
£Tl  ary  laver  on  a  surface  encounters  the  blunt  nose  of  an 
lirfoil-shaped  body  or  strut  protruding  from  that  surtace.  when 
he  approach  boundary  layer  encounters  the  wmg,  it  separates 
n  front  of  the  nose  and  the  separation  line  stretches  around 
he  sides  of  the  wing.  Closer  to  the  nose  a  horseshoe  vortex 
:orms  that  stretches  around  the  body.  A  simplified  schematic 
Df  che  horseshoe  vortex  system  is  illustrated  in  rig.  U  me 
junction  vortex  system  is  of  interest  in  many  practical  engi¬ 
neering  applications.  This  flow  pattern  is  present  at  aircraft 
wmg  roots,  ship  and  submarine  appendage/ hull  junctions,  m 
gas  turbines  at  the  blade/hub  junction,  and  m  flows  with  in- 
iecrion  normal  to  che  mean  flow. 

The  flowfield  in  che  wing/body  junction  is  dominated  by 
large-scale,  aperiodic  motions.  Devenport  and  bimpson 
measured  velocities  in  the  plane  of  symmetry  upstream  or  a 
wing/ body  junction  and  observed  bimodal  (douole  penned) 
velocity  probability  histograms.  Further  research  showed  chat 
the  flow  at  a  given  point  in  this  zone  switched  apenodi<ally 
from  one  mode  to  the  other.  This  douole-peaked  structure 
has  also  been  found  in  histograms  of  pressure  fluctuation 
measurements--’  made  on  che  endwall  upstream  or  a  wmj 
body  junction.  Hydrogen  bubble  flow  visualizations  by  Rim 
et  al./  revealed  that  the  aperiodic  motions  were  cause 
stretching  of  the  horseshoe  vortex  about  the  wing  an  ■  . 
interaction  of  the  horseshoe  vortex  with  corotat.ng  *pam»n 
vortices  that  form  between  che  horseshoe  vortex  and  ‘ 
of  separation.  This  aperiodic  phenomenon  is 
the  observed  high  turbulence  intensities  and  high  surface  pres¬ 
sure  fluctuations  and  may  be  responsible  tor  hign  neat  trans  . 
rates.  This  is  of  great  importance  to  gas  curDine  heat  .ro  • 
wnere  hub  and  casing  boundarv  layers  produce  ,ar3 
transfer  rates  at  the  functions  with  blades. 

Few  studies  of  che  effects  of  :ne  horseshoe  vortex 
.vail  heat  flux  have  been  made.  Time-averaged  heat  t  ■ 
mremencs  n  turome  cascades  nave  -een  made  ■-> 


LEWIS.  SIMPSON.  AND  UlLLtK  bUKhALc  nc.-A  <  ru- 


researchers.  Blairs  and  Dunn  and  Stoddard*  measured  end- 
wai!  heat  flux  rates  between  blades  in  a  turbine  cascade  and 
found  that  the  horseshoe  vortex  increased  heat  flux  by  a  factor 
of  3  near  the  leading  edge  of  the  turbine/endwall  junction. 

The  effects  of  the  approach  boundary-layer  thickness  on 
heat  flux  to  the  endwall  in  turbine  cascades  have  been  inves¬ 
tigated  with  inconclusive  results.  Grazianni  et  ai.7  found  that 
minimum  and  maximum  Stanton  number  levels  were  inde¬ 
pendent  of  the  approach  boundary-layer  thickness.  Grazianni 
et  al.  also  found  that  the  area  of  influence  of  the  horseshoe 
vortex  was  increased  with  increasing  approach  boundary-layer 
thickness,  in  contrast.  Georgiou  et  ai/  found  that  the  effect 
of  reducing  the  inlet  boundary-layer  thickness  was  to  uni¬ 
formly  increase  the  local  heat  flux  rate,  but  that  the  shape  of 
the  iso-heat-transfer  contours  were  unaffected  by  boundary 
laver  thickness.  Georgiou  et  al.  also  found  changes  in  tree- 
stream  turbulence  levels  to  have  only  a  minimal  effect  on  heat 
flux  rates  in  the  horseshoe  vortex. 

Gaugler  and  Russell*  compared  their  measured  turbine/ 
endwall  heat  flux  rates  with  flow  visualizations  of  Hylton  et 
al.m  The  only  obvious  correlation  between  the  horseshoe  vor¬ 
tex  and  endwall  heat  flux  was  found  near  the  vane  leading 
edge  where  a  local  peak  in  heat  flux  occurred.  The  three- 
dimensional  separation  line  was  found  to  not  correlate  with 
any  endwall  heat  flux  features.  Hippensteele  and  Russell11 
measured  turbine/endwall  heat  flux  rates  using  high  spatial 
resolution  liquid  crystal  sheets  for  the  same  turbine  cascade 
used  by  Gaugler  and  Russell  and  obtained  similar  results. 

Only  recently  have  time-resolved  heat  flux  measurements 
been  made  in  the  junction  vortex.  A  high-frequency-response 
heat  flux  gauge  called  the  Heat  Flux  Microsensor  was  recently 
used  by  Swisher  et  al.i:  to  measure  time-resolved  endwall 
heat  flux  aiong  the  stagnation  streamline  upstream  of  a 
streamlined  cylinder.  The  rms  of  the  heat  flux  unsteadiness 
was  found  to  be  as  high  as  30%  of  the  mean  heat  flux  in  the 
horseshoe  vortex. 

All  of  the  previous  studies  have  shown  similar  levels  and 
patterns  of  wing/endwall  heat  flux.  High  levels  of  heat  flux 
were  found  in  the  horseshoe  vortex  with  peak  levels  occurring 
near  the  wing.  Freestream  turbulence  and  approach  bound- 
arv-laver  thickness  have  been  found  to  have  little  effect  on 
heat  flux  to  the  endwail.  but  these  effects  are  still  poorly 
understood.  The  relationship  between  the  fluid  dynamics  and 
the  heat  flux  has  not  been  established. 

The  objective  of  this  article  is  to  relate  the  turbuient  heat 
flux  to  the  fluid  dvnamics  in  a  horseshoe  vortex.  No  attempt 


is  made  here  to  predict  the  fluid  dynamics  because  this  has 
been  the  subject  of  many  previous  papers.  It  will  therefore 
be  assumed  that  the  fluid  dynamics  (location  and  strength  of 
the  vortex)  are  known  and  the  heat  flux  will  be  related  to 
these  parameters.  The  primary  model  chosen  for  this  work 
was  the  3:2  semielliptic  nosed/NACA  0020  tailed  airfoil  for 
which  Devenport  and  Simpson1  have  made  detailed  3-com¬ 
ponent  laser  Doppler  velocimetrv  (LDV)  measurements. 

Here  we  report  detailed  time-  and  spatially-resolved  heat 
flux  measurements  that  are  necessary  to  understand  the  ef¬ 
fects  of  the  horseshoe  vortex  on  the  heat  flux.  We  examine 
the  effects  of  turbulence  on  heat  flux  and  determine  if  the 
bi modal  velocitv  phenomenon  is  responsible  for  the  augmen¬ 
tation  in  the  surface  heat  flux.  Results  for  several  different 
wing  shapes  also  are  presented  to  examine  the  effects  of  wing 
geometry. 

Experimental  Apparatus  and  Procedure 
Test  Facility 

The  wind  tunnel  is  an  open-circuit  type  and  is  powered  by 
a  centrifugal  blower.  Air  from  the  blower  is  supplied  to  the 
test  section  after  passing  through  a  fixed-setting  damper,  a 
plenum,  a  section  of  honeycomb,  seven  screens  which  are 
used  to  remove  much  of  the  turbulence  intensity,  and  a  4:1 
contraction  nozzle  to  further  reduce  the  turbulence  and  to 
accelerate  the  flow  to  test  speed.  The  potential  core  of  the 
flow  entering  the  test  section  is  uniform  to  within  0.5%  in  the 
spanwise  direction  and  1%  in  the  vertical  direction  with  a 
turbulence  intensity  of  0.1%  at  27  m/s.1 

Figure  2  is  a  side  view  of  the  6-m-iong  and  0.91-m-wide 
test  section.  The  upper  wall  is  made  from  Plexiglas^  rein¬ 
forced  with  aluminum  channel.  The  glass  side  walls  are  line 
internally  with  removable  6.4-mm-thick  Plexiglas  sheets.  The 
lower  wall  is  made  from  19-mm-thick  fin-form  plywood,  ow 
entering  the  test  section  is  subjected  to  a  further  1.0:1  con¬ 
traction'  produced  by  the  shape  of  the  upper  wail.  A  throat 
is  reached  1.63  m  downstream  of  the  entrance  where  the 
section  is  254  mm  in  height.  Downstream  of  the  throat  the 
upper  wall  is  almost  parallel  to  the  flat  lower  wail,  diverging 
gradually  from  it  with  distance  downstream  to  account  tor 
boundarv-laver  growth. 

In  order  to  create  a  new  fluid  boundary  laver  at  the  begin- 
nine  of  the  heat  flux  surface,  a  heated  false  floor  was  placed 
in  the  tunnel  and  a  suction  slot  was  oDened  in  the  tunnel  floor 
as  shown  in  Fie.  2.  The  uostream  bounoarv  layer  iormed  on 
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Modified  NACA  0020  Streamlined  Cylinder  NACA  0015 


Fig.  1  Side  view  of  the  boundary-layer  tunnel  test  section  with  heated  false  floor  and  wing  shapes. 


the  tunnel  floor  is  sucked  ouc  through  the  slot  and  a  new 
boundary  layer  develops  from  the  rounded  leading  edge  of 
the  false  floor.  The  thin  boundary  layer  is  tripped  by  a  12.7- 
mm-wide  strip  of  120-grade  sandpaper,  located  with  its  lead¬ 
ing  edge  70  mm  downstream  of  the  false  floor  leading  edge. 
This  is  the  same  arrangement  used  by  Oilmen  and  Simpson.1 
The  forward  edge  of  the  suction  slot  is  located  145  mm  down¬ 
stream  of  the  throat  of  the  test  section.  The  false  floor  rests 
with  its  upper  surface  35  mm  above  the  tunnel  floor  and  its 
leading  edge  116  mm  downstream  of  the  throat.  The  false 
floor  was  made  of  3  sections  of  1.59-cm-thick  aluminum  plate. 
A  craiiing-edge  flap  set  at  an  angle  of  7  deg  was  used  to 
produce  a  constant  static  pressure  on  the  tunnel  floor  through¬ 
out  the  heated  test  section. 

The  wing  was  mounted  in  the  test  section  at  zero  incidence 
and  sweep  with  its  leading  edge  1.17  m  downstream  of  the 
leading  edge  of  the  false  floor.  Since  the  measurements  were 
made  only  with  the  wing  at  a  zero  angle  of  attack,  it  was 
possible  to  reduce  blockage-induced  pressure  gradients  by 
contouring  the  side  walls  of  the  wmd  tunnel  to  approximately 
follow  the  streamlines  produced  by  the  3:2  elliptic  nose/NACA 
0020  tailed  wing  in  unbounded  potential  flow.  This  was  ac¬ 
complished  by  removing  the  6-,4-mm-thick  Plexiglas  sheets 
that  line  the  tunnel  sidewalls  from  the  region  surrounding  the 
wing.  This  effectively  increased  the  width  of  the  test  section 
by  12.7  mm  from  a  location  330  mm  upstream  of  the  wing 
leading  edge  to  another  203  mm  downstream  of  its  trailing 
edge.  The  abrupt  comers  at  the  edges  of  the  liner  were  faired 
over  using  adhesive  tape.  This  same  arrangement  was  used 
for  all  wing  shapes. 

Streamwise  velocity  measurements  were  made  using  a  hoc- 
wire  anemometer  to  determine  the  cold  wall  approach  bound¬ 
ary-layer  conditions.  The  approach  boundary-layer  velocity 
profile  was  measured  21.58  cm  in  front  of  the  leading  edge 
of  the  3:2  elliptic  nose/NACA  0020  tailed  wing  with  the  model 
in  place.  The  nominal  freestream  velocity  was  32.34  m/s  with 
a  Reynolds  number  based  on  momentum  thickness  of  3730. 
The  momentum  thickness  was  L.9  mm  and  the  displacement 
thickness  was  2.6  mm.  The  friction  velocity  was  estimated  to 
be  1.32  m/s. 

The  aluminum  false  floor  was  heated  and  held  at  a  constant 
and  uniform  temperature.  Silicone-rubber-insulated  electric 
resistance  heaters  were  held  against  the  bottom  of  the  false 
floor  by  l.9-cm-(hick  styrofoam  sheet  which  was  secured  by 
bolts.  T-type  thermocouples  were  soldered  into  3.  i8-mm-diam 
brass  tube  and  press  fit  into  the  aluminum  floor  for  feedback 
control.  Power  to  the  heaters  was  controlled  by  5.  Surotherm 
510.  3-mooe  process  controllers,  each  connected  to  a  £u- 
rotherm  531  SCR  power  suoply.  Because  of  the  high  thermal 
conductivity  of  aluminum  and  the  thickness  of  the  raise  floor. 


the  surface  was  nearly  isothermal.  The  automatic  controllers 
held  the  surface  temperature  constant  and  uniform  to  within 
-0.5°C.  The  location  of  a  ••virtual  origin”  was  calculated  from 
the  momentum  thickness  using  the  correlation  for  turbulent 
boundary  layers  given  by  Kays  and  Crawford1-*: 

(dfx)  =  0.03 6Re;"z  (0 

The  virtual  origin  was  i  cm  downstream  of  the  leading  edge 
of  the  false  floor.  The  unheaced  starring  length  was  approx¬ 
imately  5  cm.  and  was  therefore  negligible. 

The  heat  flux  probe  was  mounted  in  a  cam  system  installed 
in  the  false  floor  in  front  of  the  wing.  The  cams  were  made 
of  1.59-cm-thick  aluminum  plate  and  were  heated  from  be¬ 
neath  by  electric  resistance  heaters.  Rotating  one  or  bothj:ams 
allowed  the  probe  to  be  positioned  throughout  a  L2.7-cm- 
radius  region  in  front  of  the  wing.  The  cam  system  is  illustrated 
in  Fig.  7  of  Lewis  et  al.u 

Wing  Shapes 

The  wing  shapes  used  were  a  modified  NACA  0020  tailed 
wing,  a  streamlined  cylinder  shape,  and  an  NACA  0015.  The 
modified  NACA  0020  wing  consists  of  a  3:2  elliptical  nose 
(with  its  major  axis  aligned  with  the  chord)  and  an  NACA 
0020  tail  joined  at  the  maximum  thickness.  These  wing  shapes 
were  chosen  because  of  the  bimodai  surface  pressure  fluc¬ 
tuation  phenomenon  observed  by  Oilmen  and  Simpson-  for 
these  wing  shapes.  The  height  of  each  model  was  2-.9  cm. 
which  spanned  the  test  section.  The  models  were  each  made 
of  aluminum  and  were  heaced  with  silicone-rubber-insulated 
electric  resistance  heaters  glued  inside  the  models.  Each  model 
was  positioned  in  the  same  location  in  the  tunnel  with  its 
leading  edge  1.17  m  from  the  leading  edge  of  the  false  floor. 
The  flow  over  each  model  was  tripped  with  a  strip  of  120- 
grade  sandpaper  to  prevent  unsteadiness  due  to  natural  flow 
Transition.  The  length,  maximum  thickness,  trip  width,  and 
tnp  location  of  each  wmg  shape  are  given  in  Table  l.  Results 
from  oil  flow  visualizations  made  by  dlcmen  and  Simpson  on 
the  false  floor  in  front  of  each  wing  are  also  given  in  Table 
l.  These  oil  visualizations  revealed  the  location  of  primary 
separation  of  the  approach  boundary  laver  and  the  lo¬ 
cation  of  a  second  line,  closer  to  the  model  which  Devenport 
and  Simpson1  have  shown  to  mark  the  location  of  minimum 
local  streamwise  shear  c,(>u  measured  from  the  wing  leading 
edge  along  the  line  of  symmetry  in  front  of  each  wing.  These 
lines  are  illustrated  in  Fig.  3. 

Heat  Flux  Measurement 

Heat  flux  measurements  were  made  using  both  a  hign-(re- 
quencv-resDonse,  low  outDut  gauge  ana  a  nign  sensinvitv  gauge. 


Lfc.  W 1  /Ai’iw  rxi  Li  c  rv  ; 


■; . . . - 


Table  I  Characteristic  dimensions  of  the  wing  models  and  results  from  oil  flow  visualization 


Model  name 

Chord 

length. 

cm 

Maximum 

thickness, 

cm 

Trip 

location/chord. 

circumferential 

Trip  length/ 
chord 

W 

cm 

Modified 

NACA  0020 
Streamlined 
cylinder 

NACA  0015 

30.5 

29.8 

60.8 

7.17 

12,70 

9.21 

0.1406 

0.2542 

0.1152 

0.0208 

0.0374 

0.0376 

3.23 

6.59 

2.71 

0.60 

0.53 

0.53 

•‘From  RcL  2. 


x/ *seo 


Fig.  3  Contours  of  constant  time-mean  Stanton  number  x  103  for 
the  modified  NACA  0020  wing;  — .  line  of  separation;  +  +  +  +  +■ 
line  of  low  shear:  •,  locations  of  heat  flux  measurements.  Line  1  is 
along  the  line  of  symmetry,  with  lines  2-7  proceeding  clockwise  around 
the  wing. 

The  hich-frequencv-response  gause  was  a  Heat  Flux  Micro- 
sensor' “manufactured  by  Vatell  Corp-  Similar  gauges  have 
been  used  in  several  other  studies  of  turbulent  heat  flux. 

The  Heat  Flux  Microsensor  is  a  layered-type  gauge  that  op¬ 
erates  bv  relating  the  heat  flux  to  the  measured  temperature 
drop  across  a  thin  thermal  resistance  layer.  A  detailed  de¬ 
scription  of  the  Microsensor  construction  can  be  found  m 
Hager  et  al.,f>  The  Microsensor  has  a  flat  frequency  response 
up  to  50  kHz'7  and  outputs  a  voltage  directly  proportional  to 
the  heat  flux.  The  gauge  was  calibrated  by  measuring  the  heat 
flux  rate  at  a  three-dimensional  stagnation  point  where  the 
heat  flux  coefficient  had  been  previously  determined.  The 
gauge  sensitivity  was  found  to  be  66  -  10  /xV/(W/cm-).‘*  The 
thin^  construction  of  the  gauge  (<2  ^m)  is  physically  unob¬ 
trusive  in  the  flow  and  causes  minimal  disruption  of  the  ther¬ 
mal  boundary  laver  due  to  a  step  change  in  temperature.  The 
sensing  area  of  the  gauge  is  3  x  3  mm.  Two  0.076-mm-diam 
tvpe-T  thermocouples  were  mounted  in  the  surface  of  .the 
gauge  and  the  floor  to  match  the  gauge  surface  temperature 
with  the  temperature  of  the  heated  floor. 

The  output  voltage  of  the  Microsensor  was  amplified  by  a 
differential  dc-powered  amplifier  designed  and  built  by  Vatell 
Corp.  The  signal  was  amplified  by  a  gain  of  1000  and  con¬ 
nected  to  a  Hewlett  Packard  3562A  Dynamic  Signal  Analyzer. 
The  signal  was  sampled  at  a  frequency  of  10.240  Hz.  One 
hundred  loads  of  2048  samples  were  taken  over  a  period  of 
approximately  1  min  and  recorded  on  floppy  disk. 

Due  to  the  low  output  voltage  of  the  Microsensor,  a  second, 
high-sensitivity  Schmidt-Boelter  gauge  manufactured  by 


Medtherm  Corp.  (model  8-2-.625-36-20893T)  was  used  to 
measure  the  mean  time-averaged  heat  flux.  The  Schmidt- 
Boelter  gauge  has  an  output  voltage  directly  proportional  to 
the  heat  flux  with  a  sensitivity  of  3.9  z:  0.2  mV/(W/cm-)  and 
a  response  time  of  approximately  250  ms.  The  diameter  of 
the  sensing  area  of  the  gauge  is  3.18  mm.  The  output  voltage 
from  the  Schmidt-Boelter  gauge  was  amplified  1000  times  and 
filtered  by  a  300-Hz  low-pass  filter.  The  signal  was  sampled 
at  a  frequencv  of  600  Hz  and  recorded  by  an  IBM  PC  AT 
using  a  12-bit  Data  Translation  DT2801A  A/D  convener.  The 
temperature  of  the  gauge  was  measured  by  a  T-type  ther¬ 
mocouple  mounted  inside  the  gauge  by  the  manufacturer. 

Heat  transfer  coefficients  were  calculated  from  the  mea¬ 
sured  gauge  output  by 

<2) 

The  convection  heat  transfer  coefficients  were  converted  to 
nondimensional  Stanton  numbers,  where  the  freestream  ve¬ 
locity,  C„,  and  p  were  evaluated  at  the  film  temperature.  The 
emissivity  of  the  Microsensor  was  approximately  0.  The 
emissivitv  of  the  Schmidt-Boelter  gauge  was  0.97. 

Experimental  Uncertainties 

Uncertainties  in  the  Stanton  numbers  were  estimated  at 
20:1  odds  using  the  methods  described  in  Kline  and  Mc- 
Clintock.1''  The'major  contributor  to  the  uncertainty  m  Stan¬ 
ton  number  was  the  uncertainty  in  the  static  sensitivity  of  the 
heat  flux  °auees.  The  uncertainties  in  Stanton  number  were 
estimated  to  'be  5.2  and  15.2%  using  the  Schmidt-Boelter 
gauge  and  the  Heat  Flux  Microsensor,  respectively.  Uncer¬ 
tainty  in  the  sensor  coordinates,  x,  z  was  0.5  mm. 

Experimental  Results  and  Discussion 

Experimental  results  are  presented  using  the  coordinate 
svstem  shown  in  Fig.  1.  Lengths  are  normalized  using  the 
distance  from  the  model  leading  edge  to  the  location  of  the 
primarv  separation  line  as  measured  from  oil-flow  visualiza¬ 
tions  and  presented  in  Table  1.  Note  that  this  distance  is 
different  for  each  model.  For  the  3:2  elliptic  nose/NACA  0020 
tailed  model,  measurements  were  taken  at  the  60  locations 
along  the  seven  lines  shown  in  Fig.  3.  The  lines  are  numbered 
clockwise  from  the  line  of  symmetry.  Lines  1,  4.  o.  and  / 
correspond  to  the  same  locations  where  Devenport  and 
Simpson1  made  3-component  velocity  profile  measurements. 
Measurements  forthe  other  two  models  were  taken  only  along 
the  centerline  upstream  of  the  model  at  the  same  locations 
where  Olpmen  and  Simpson-  measured  pressure  fluctuations. 

Mean  Heat  Flux 

Mean  heat  flux  measurements  presented  here  were  ob- 
tamed  from  the  Schmidt-Boelter  gauge.  Contours  of 
Stanton  number  are  shown  in  Fig.  3  for  the  modified  NACA 
0020  model.  The  Stanton  numbers,  nondimenstonalized  on 
freestream  conditions,  represent  lines  of  constant  heat  trans¬ 
fer  coefficient.  The  region  of  the  flow  close  to  tne  wins  is  an 
area  of  high  heat  flux  created  bv  the  presence  of  the  horseshoe 
vortex.  This  high  heat  flux  is  the  result  of  both  hignlv  turbu- 


lent  mixing  created  by  the  vortex  and  cold  freestream  air 
transported  by  the  vortex  down  the  front  of  the  wing  perpen¬ 
dicular  to  the  endwail.  The  locations  ot'  maximum  levels  of 
heat  flux  were  dose  to  the  wing,  at  the  wing/'endwall  junc¬ 
tion.  Upstream  of  the  wing,  the  Stanton  numbers  approach 
the  two-dimensional,  zero  pressure  gradient  boundary-layer 
Stanton  number  Sc  =  0.00184.  calculated  using  ScPr*A  = 
0.0287/?^  ,r 0':  and  the  estimated  location  of  the  virtual  origin. 
The  unheated  starting  length  was  negligible. 

Distributions  of  mean  surface  heat  tlux  on  the  centerline 
upstream  of  each  model  are  shown  in  Fig.  4.  From  examining 
Fig.  4,  it  is  apparent  that  the  heat  flux  enhancement  begins 
close  to  the  separation  line  (x/xsep  =  -1.0)  for  each  model 
tested  and  that  the  heat  flux  in  this  region  is  produced  by  the 
secondary  flow.  Downstream  of  separation,  the  heat  flux  in¬ 
creases  rapidly  as  the  wing  is  approached,  reaching  a  level 
almost  100%  higher  than  in  the  approach  boundary  layer. 
This  initial  rise  is  followed  by  a  plateau  in  the  heat  flux  dis¬ 
tributions.  The  heat  flux  levels  attained  here  are  approxi¬ 
mately  the  same  for  the  three  models  tested.  For  the  modified 
NACA  0020  and  the  NACA  0015  models,  this  plateau  appears 
to  be  centered  around  the  line  of  low  shear,  atx/x^,  =  -  0.60 
and  xix ^  =  -0.53.  respectively.  An  inflection  point  in  the 
plot  of  Stanton  number  vs  x/x^p  occurs  just  downstream  of 
the  line  of  low  shear,  approximately  below  the  vortex  center, 
followed  by  a  second  rapid  increase  in  heat  flux.  The  maxi¬ 
mum  value  of  Stanton  number  was  found  at  the  closest  mea¬ 
surement  location  to  the  wing,  where  the  maximum  Stanton 
number  measured  was  0.0054,  a  200%  increase  over  the  value 
of  0.0018  measured  in  the  approach  boundary  layer.  For  the 
tapered  cylinder  model,  the  plateau  in  heat  flux  extends  into 
the  wing/endwall  junction  with  no  second  rise  in  heat  flux. 

Heat  Flux  Fluctuations 

Fluctuations  in  endwail  surface  heat  flux  presented  here 
were  measured  using  the  Heat  Flux  Microsensor.  The  fluc¬ 
tuating  component  of  heat  flux  measured  along  the  centerline 
in  front  of  each  model  is  shown  in  Fig.  5.  The  fluctuations  in 
heat  flux  increase  rapidly  after  the  separation  line  and  reach 
a  maximum  value  near  the  line  of  low  shear.  Fluctuations  in 
heat  flux  decrease  rapidly  as  the  wing  is  approached.  The 
fluctuating  heat  flux  is  a  significant  part  of  the  total  heat  flux, 
reaching  values  as  high  as  25%  of  the  mean  heat  flux  near 
the  line  of  low  shear.  The  line  of  low  shear  has  also  been 
found  ta  be  the  location  of  maximum  level  of  pressure 
fluctuations--*  and  turbulence  normal  stresses.1 

Probability  Density  Functions 

Distributions  of  the  probability  density  function  (pdf)  of 
the  heat  flux  fluctuations  measured  along  the  centerline  in 
front  of  the  modified  NACA  0020  airfoil  are  shown  in  Fig. 
6.  The  pdfs  of  heat  flux  fluctuations  do  not  show  the  bimodal 
or  “double-peaked”  structure  that  have  been  found  in  the 


Fig.  4  Time-mean  Stanton  number  distribution  along  the  stagnation 
streamline  tor  each  model. 
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Fig.  5  Fluctuations  in  Stanton  number  measured  along  the  stagnation 
streamline  for  each  model. 
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Fig.  6*  Probability  density  functions  of  heat  flux  fluctuations  mea¬ 
sured  on  the  stagnation  streamline  for  the  modified  NACA  0020  wing. 


velocity  pdfs  and  in  pdfs  of  pressure  fluctuations  measured 
on  the  endwail  in  front  of  the  wing. 

Although  the  heat  flux  pdfs  are  not  bimodal.  the  effects  of 
large  scale  unsteadiness  are  seen  in  the  distortion  of  the  pdfs. 
Distributions  of  the  skewness  and  flatness  factors  of  the  heat 
flux  pdfs  are  shown  for  each  model  in  Fig.  7  and  Fig.  S, 
respectively.  For  reference,  a  Gaussian  pdf  has  a  skewness 
equal  to  0  and  a  flatness  factor  equal  to  3.  Maximum  values 
of  skewness  and  flatness  of  the  heat  flux  pdfs  were  round  in 
the  vicinity  of  separation  and  are  caused  by  the  highly  inter¬ 
mittent  forward  and  reverse  flow  which  exists  in  this  region. 
Table  2  presents  the  locations  and  values  of  the  maximum 
mean  Stanton  number,  rms  Stanton  number,  skewness,  ana 


Table  2  Location  and  value  of  maximum  and  minimum  5L  St skewness,  and  flatness 


Modified 

NACA  0020 

Streamlined 

cylinder 

NACA  0015 

Value 

x/x^r 

Value 

xixap 

Value 

n 

Maximum  St  x  10' 

5.37 

-0.118 

3.51 

-0.131 

4.77 

-0.078 

Maximum  St m,s  x  10' 

0.754 

-0.636 

0.833 

-0.498 

0.686 

-0.524 

Maximum  skewness 

1.08 

-0.909 

1.23 

-0.876 

0.96 

-0.898 

Maximum  flatness 

4.82 

-i.or 

5.61 

-0.946 

4.53 

-0.973 

X/ X^op 


Fig,  7  Skewness  factor  of  the  heat  flux  pdfs  for  each  model. 


x/ Xa«p 


Fig.  8  Flatness  factor  of  the  heat  flux  pdfs  for  each  model. 

flatness  for  each  wing.  Upstream  of  the  wing,  in  the  approach 
boundary  layer,  the  flatness  approaches  the  Gaussian  value 
of  3  and  the  skewness  factor  is  approximately  0.35.  The  pos¬ 
itive  skewness  in  the  approach  boundary  layer  may  be  caused 
by  the  '‘sweeps''  of  high  momentum,  cold  fluid  from  the  outer 
region  which  occurs  in  the  inner  wall  layer. 

Power  Spectra 

One-sided  autospectrai  density  functions  (power  spectra) 
of  the  heat  flux  fluctuations  were  computed.  Figure  9  presents 
a  comparison  of  nondimensionai  spectra  of  the  heat  flux  fluc¬ 
tuations  measured  in  the  approach  boundary  layer  ( xix - 
-1.4-65),  in  the  region  between  the  separation  line  and  the 
vortex  center  {x/x^  =  -0.659),  and  near  the  wing/endwai! 
junction  (jc/jr^p  =  -0.189)  for  the  modified  NACA  0020 
wing.  Frequency  has  been  normalized  by  t/Urcf.  which  is  the 
time  it  takes  the  freestream  flow  to  travel  one  wing  thickness. 
The  spectra  have  been  normalized  by  ( Urc{It)(pUrefCnlT)z  so 
that  the  area  under  the  spectrum  is  equal  to  the  variance  of 
the  local  Stanton  number.  The  spectra  have  been  nondimen- 
sionalized  in  this  way  to  allow  direct  comparison  of  spectral 
levels  measured  at  different  streamwise  locations.  Lines  in¬ 
dicating  the  -  1  and  -5/3  power-law  slopes  for  velocity  fluc¬ 
tuation  spectra*11  and  a  -8/3  power-law  spectra  are  included 
in  Fig.  9  for  comparison. 


10-2  10-1  10°  101 


ft/Urf 

Fig.  9  Comparison  of  heat  flux  spectra  measured  in  the  approach 
boundary  layer  (x/x^  =  - 1.465),  in  the  region  between  the  separation 
line  and  the  line  of  low  shear  {xix^  =  -0.659),  and  near  the  junction 
ixlx^  =  -0.189)  for  the  modified  NACA  0020  wing. 

In  general,  the  spectra  have  a  region  of  almost  constant 
spectral  level  at  low  frequencies  followed  by  a  rolloff  with  a 
-8/3  slope.  While  no  spectra  of  heat  flux  fluctuations  were 
available  in  the  literature,  the  -8/3  power  law  has  been  found 
in  the  spectra  for  temperature  fluctuations  measured  in  fully 
developed  pipe  flow. 21  Note  that  the  results  have  not  been 
corrected  for  the  effects  of  the  finite  size  of  the  heat  flux 
sensor  and  that  the  spectra  at  high  frequencies  may  be  atten¬ 
uated.  Upstream  of  separation,  in  the  approach  boundary 
layer,  the  region  of  near  constant  spectral  level  extends  up  to 
a  dimensionless  frequency  of  ft!Unt  =  0.3,  followed  by  the 
-8/3  rolloff.  Downstream  of  separation,  between  the  line  of 
separation  and  the  vortex  center,  there  is  a  dramatic  increase 
in  spectral  density,  particularly  at  low  frequencies.  Further 
downstream,  between  the  vortex  center  and  the  wing/endwall 
junction,  there  is  an  increase  in  spectral  density  at  high  fre¬ 
quencies. 

The  high  levels  of  spectral  density  found  in  the  region  be¬ 
tween  the  separation  line  and  the  vortex  center  are  associated 
with  the  high  heat  flux  fluctuations  shown  in  Fig.  5.  The  high 
spectral  levels  at  low  frequencies  suggest  that  the  heat  flux 
fluctuations  in  this  region  are  created  by  large-scale,  low- 
frequency,  aperiodic,  coherent  motions.  Spectra  of  the  nor¬ 
mal  and  streamwise  velocity  components  measured  by  a  laser- 
Doppler  velocimeter  in  this  region show  a  similarly  dramatic 
increase  in  spectral  density  at  low  frequencies. 

Near  the  wing/endwail  junction,  spectral  energy  exists  at 
higher  frequencies  than  elsewhere  in  the  flowfield.  At  low 
frequencies,  the  spectral  density  levels  are  the  same  as  in  the 
approach  boundary  laver.  However,  the  region  of  constant 
SDectrai  density  extends  further,  to  a  dimensionless  frequency 
value  of  ft/Unl  =  i.O.  This  high-frequency  energy  could  be 
the  effect  of  a  small  counter-rotating  vortex  trapped  in  the 
junction,  or  the  effect  of  small  structures  from  the  outer 
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boundary-layer  region,  transported  by  the  main  vortex,  down 
the  wing  and  into  the  junction. 

Fromthe  results  presented  above  and  the  available  velocity 
data,  a  description  of  the  augmentation  in  heat  flux  across 
the  horseshoe  vortex  can  be  summarized  by  breaking  the 
flowfield  into  three  regions.  In  the  first  region,  upstream  of 
the  line  of  separation,  the  heat  flux  is  as  would  be  expected 
for  a  two-dimensional,  adverse  pressure  gradient  boundary 

layer.  .  ,  . 

The  second  region,  between  the  line  of  separation  ana  the 
line  of  low  shear,  is  dominated  by  the  large-scale  aperiodic 
motions  of  the  horseshoe  vortex.  This  region  is  characterized 
by  highly  intermittent  forward  and  reverse  flow.  Turbulence 
energy  production  and  turbulent  stresses  an  order  of  mag¬ 
nitude  higher  than  upstream  were  found  in  this  region  al¬ 
though  surface  shear  stresses  are  small.1  Mean  heat  flux  levels 
are  increased  over  upstream  values  by  a  factor  of  2,  and 
fluctuations  of  the  heat  flux  are  increased  by  a  factor  of  3. 
The  augmentation  in  heat  transfer  in  this  region  is  due  to 
turbulent  mixing  created  by  the  unsteadiness  of  the  recircu¬ 
lating  flow. 

Downstream  of  the  line  of  low  shear,  in  the  third  region, 
the  heat  flux  increases  rapidly  and  reaches  a  maximum  value 
in  the  immediate  vicinity  of  the  wing/endwall  junction  while 
turbulence  energy  and  fluctuations  in  heat  flux  decrease  rap¬ 
idly  as  the  wing  is  approached.  Skin  friction  values  in  this 
region  are  very  large.12  The  velocity  field  in  this  region  is 
similar  to  an  impinging  jet.  The  fluid  transported  toward  the 
endwall  on  the  downwash  side  of  the  vortex  impinges  on  the 
wing/endwall  junction  and  is  deflected  radially  in  the  spanwise 
direction  and  back  upstream  in  a  near  wail  jet  of  fluid.  The 
heat  flux  in  this  region  is  determined  by  the  effects  of  the 
vortex  on  the  mean  flow.  The  increase  in  heat  flux  in  this 
third  region  is  due  to  a  thinning  of  the  thermal  boundary  layer 
as  cold  fluid  from  the  outer  boundary  layer  or  freestream  is 
transported  by  the  vortex  to  the  near  wail  region. 

This  may  explain  the  differences  in  mean  heat  flux  in  the 
near-wing  region  for  each  wing  shown  in  Fig.  4.  The  time- 
mean  heat  flux  in  this  near-wing  region  appears  to  be  con¬ 
trolled  by  the  size  and  location  of  the  primary  vortex.  The 
size  of  the  primary  vortex  is  proportional  to  the  location  of 
the  line  of  separation  which  was  3.23.  2.71.  and  6.59  cm 
for  the  modified  NACA  0020,  NACA  0015.  and  streamlined 
cylinder  models*  respectively.  For  the  streamlined  cylinder, 
the  vortex  is  located  much  farther  upstream,  away  from  the 
wing  than  for  the  other  two  models.  This  may  lead  to  a  lower 
downwash  velocity  in  the  near- wing  region  for  the  streamlined 
cylinder  model.  Unfortunately,  no  detailed  velocity  data  is 
available  for  this  model  to  verify  this  conclusion. 

Evaluation  of  SC  Correlation 

Since  the  heat  flux  rate  outside  of  the  line  of  low  shear 
(in  the  first  and  second  regions  defined  above)  is  believed  to 
be  primarily  due  to  turbulence  generated  by  the  large-scale 
unsteadiness  of  the  horseshoe  vortex,  one  would  expect  a 
correlation  between  the  turbulent  stresses  and  the  time-av¬ 
eraged  heat  tlux.  For  two-dimensional  boundary-layer  flows. 
Maciejewski  and  Moffat23  relate  the  local  heat  tlux  to  the 
streamwise  component  of  the  turbulent  normal  stress  through 
the  use  of  a  Stanton  number  based  on  as 

SC  ^  /t/(pC,urTO.mjJ  (3) 

where  h  is  the  local  convection  coefficient  and  ^  is  the 
maximum  value  of  urms  found  in  the  boundary  layer  or  tree- 
stream  at  the  location  where  ft  is  measured.  Mote  that  St'  is 
not  to  be-  confused  with  which  is  the  nondimenstonal 

rms  value  of  the  fluctuations  in  heat  tlux.  The  St'  correlation 
is  of  limited  practical  value  since  it  requires  a  detailed  knowl¬ 
edge  of  the  turbulence  structure  chat  is  usually  not  available. 
The  authors  use  the  SC  correlation  only  L)  to  demonstrate 
that  the  heat  tlux  in  this  region  is  due  to  the  turbulence  rieia 


and  2)  to  demonstrate  that  the  SC  correlation  can  be  applied 
even  to  this  complex,  highly  three-dimensional,  separated  flow. 

Maciejewski  and  Moffat23  suggest  a  functional  form  of  the 
relation  between  SC  and  turbuFence  intensity  Tu  =  ^Nns.mux^ 
UrKf  for  air  (Pr  =  0.71)  as 

SC  =  0.0184  4-  0.0092  exp  -  [(Tu  -  0.1i)/0.055]2  (4) 

This  form  has  a  peak  value  at  a  turbulence  level  of  11%.  For 
turbulence  intensities  above  0.20.  SC  0.0184.  Maciejewski 
and  Moffat  applied  this  correlation  to  heat  flux  experiments 
for  various  complex  flow  situations  from  the  literature  and 
found  Eq.  (4)  to  be  valid  to  within  15%  at  (20:1)  odds  in¬ 
dependent  of  flow  geometry  or  Reynolds  number. 

To  extend  the  SC  correlation  to  three-dimensional  flows, 
the  authors  chose  to  substitute  the  maximum  value  of  tur¬ 
bulent  kinetic  energy  present  in  the  boundary  layer  A:™.™* 
in  place  of  in  Eq.  (3).  The  value  of  U was 

obtained  from  the  LDV  data  of  Devenport  and  Simpson.1 
which  were  measured  at  a  freestream  velocity  of  27  m/s  and 
an  approach  Reynolds  number  based  on  momentum  thick¬ 
ness  of  6700.  This  value  of  may  not  be  the  same 

as  the  value  of  for  the  present  data  for  which 

U^i  =  32.3  m/s  and  Re*  =  3700.  It  is  assumed*  however  that 
is  not  a  stroni  function  of  Reynolds  number  for 
this  fiow.  This  assumption  is  supported  by  the  fact  that  the 
locations  of  the  line  of  low  shear  and  the  line  of  separation 
revealed  by  the  oil  flow  visualizations  performed  at  the  lower 
Reynolds  number2  were  found  to  agree  with  the  locations  of 
these  lines  obtained  from  oil  flow  visualizations  of  Devenport 
and  Simpson.1  This  suggests  that  the  size  and  location  of  the 
primary  vortex,  and  presumably  the  turbulence  structure,  is 
the  same  for  the  flowfield  of  Devenport  and  Simpson  and  the 
flowfield  of  the  current  study. 


Fig.  10  Relation  between  beat  flux  and  in  the  region  outside 

of  the  line  of  low  shear  for  the  modified  lNACA  0020  wing. 
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Fig.  1 1  SC  correlation  for  beat  tlux  in  the  region  outside  of  the  iine 
„i  low  shear  for  the  modified  NACA  0020  wuig. 


Figure  10  shows  the  local  heat  transfer  coefficient  as  a  func¬ 
tion  of  k  .  IU",  for  data  measured  along  lines  1.  J.  6  and 
7  Onlv  data  from  reeions  I  and  II.  outside  of  the  line  of  low 
shear,  were  used.  The  heat  flux  in  this  region  vanes  by  a 
factor  of  2.5.  A  linear  relation  is  seen  to  exist  between  St  and 
/.  iu  St'  correlations  were  computed  from  the  heat 
transfer  coefficients  in  Fig.  10.  and  are  shown  as  a  function 
oft.  ju  ,  in  Fig.  11.  The  values  of  St'  have  a  mean  value 
of  0.021? and  a  standard  deviation  of  11%.  A  peak  in  St'  does 
seem  to  appear  around  a  turbulence  level  of  10%. 

Conclusions 

Time-  and  spatiallv-resolved  heat  flux  measurements  have 
been  made  on  the  e'ndwall  surface  in  a  wing/body  junction 
vortex  svstem.  Both  the  time-averaged  and  fluctuating  com¬ 
ponents  of  the  heat  flux  have  been  examined.  The  ume-av- 
eraged  heat  flux  was  found  to  have  a  maximum  value  m  the 
immediate  vicinity  of  the  wing.  This  high  level  of  heat  flux  in 
the  region  adjacent  to  the  wing  is  due  to  the  effect  of  the 
junction  vortex  on  the  mean  velocity  field.  The  position  of 
the  primary  vortex  center,  and  therefore  the  wing  shape,  was 
found  to  stronglv  affect  the  level  of  heat  flux  in  this  region. 

Further  upstream,  between  the  separation  line  and  the  line 
of  low  shear,  the  time-averaged  heat  flux  was  found  to  be 
related  to  the  high  levels  of  turbulence  stresses  produced  by 
the  bimodal  switching  of  the  velocity  field.  The  heat  flux 
measured  in  this  region  was  related  to  the  turbulent  stresses 
using  a  three-dimensional  extension  of  the  St'  correlation  of 
Macieiewski  and  Moffat. 

It  was  found  that  the  line  of  low  shear  and  the  maximum 
value  of  fluctuations  in  heat  flux  occur  at  the  same  lo¬ 
cation.  Histograms  of  the  heat  flux  fluctuations  were  highly 
distorted  in  the  region  surrounding  the  separation  line  al¬ 
though  no  bimodal  structure  was  found.  Maximum  values  of 
skewness  and  flatness  of  the  histograms  of  surface  heat  flux 
were  located  near  the  line  of  separation.  Spectra  of  the  heat 
flux  fluctuations  showed  high  levels  of  low-frequency  energy 
in  this  region. 
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